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ABSTRACT 

This  study  is  concerned  with  effects  of  diffusa  target  scintillation 
on  10.6ym  heterodyne  line  scan  systems  where  the  objective  was  to  identify 
the-'problem  areas  which  would  eventually  control  the  limiting  factors  in  the 
image  and  MTI  performance  of  line  scan  systems.  In  particular,  diffuse  tar- 
get scintillation  (speckle)  was  recognized  as  one  of  the  limiting  factors  for 
a heterodyne  receiver  design  and  this  study  is  concerned  with  the  statistics 
of  speckle  noise  and  with  speckle  reduction  techniques. 
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1 . INTRODUCTION  AND  SUMMARY 

The  following  report  is  a summary  of  the  major  topic  of  investigation 
for  the  laser  line  scanner  study  conducted  for  the  HOWLS  program  from  March 
1,  1975  to  August  1,  1975.  The  study  is  primarily  concerned  with  effects 
of  diffuse  target  scintillation  of  10.6um  heterodyne  line  scan  systems.  An 
objective  of  this  investigation  was  to  identify  the  problem  areas  which 
would  eventually  control  the  limiting  factors  in  the  image  and  MTI  perfor- 
mance of  line  scan  systems.  In  particular,  diffuse  target  scintillation 
(speckle)  was  recognized  as  one  of  the  limiting  factors  for  a heterodyne 
receiver  design  and  the  major  emphasis  of  this  report  is  concerned  with  the 
statistics  of  speckle  noise  and  with  speckle  reduction  techniques.  Speckle 
reduction  is  especially  important  for  heterodyne  systems  since,  for  single 
detector  receivers,  aperture  averaging  is  not  advantageous. 

Coherent  receiver  technology  has  been  developed  to  the  point  where  the 
advantages  of  active  coherent  imagery  and  MTI  have  been  clearly  demonstrated. 
At  10.6ym  these  systems  have  qualitatively  demonstrated  the  ability  to  ob- 
tain high  resolution  imagery,  pointing  and  tracking,  ranging,  and  doppler 
capability.  However,  quantitative  data  on  the  ultimate  performance  of 
coherent  systems  is  scarce  and  there  are  many  design  areas  which  are  in  need 
of  a quantitative  measurements  program  which  would  provide  more  information 
for  future  system  design. 

It  has  become  clear  that  some  of  the  limiting  factors  in  active  systems 
at  10.6pm  are  the  scintillations  from  diffuse  targets  (speckle),  specular 
targets  (glint)  and  atmospheric  turbulence.  The  fluctuating  signals  from 
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these  phenomena  raise  image  contrast  and  lower  the  system  S/N,  angular  re- 
solution, and  target/MTI  detection  probabilities.  In  the  report  that  fol- 
lows we  wish  to  outline  a measurement  program  that  would  obtain  i<-<»ntitative 
data  from  laboratory  measurements  which  would  enhance  future  systems  which 
are  designed  to  reduce  speckle. 

We  are  concerned  with  the  quality  of  coherent  images  obtained  from  an 
active  line  scanner  system  which  operates  at  lO.bvim.  We  wish  to  obtain  the 
maximum  quantitative  information  from  targets  and  backgrounds  of  known  re- 
flectivity and  surface  characteristics.  In  particular  we  wish  to  make  known 
the  relation  between  image  contrast  and  angular  resolution  for  coherent 
receivers  in  situations  which  the  contrast  is  dominated  by  diffuse  target 
scintillation. 

The  noise  equivalent  reflectivity  of  an  active  line  scanner  is  the 
minimum  difference  in  reflectivity  observed  and  is  limited  by  the  system 
noise  or  target  scintillation.  It  will  be  shown  that  the  background  and 
targets  diffuse  surface  qualities  can  introduce  an  associated  noise  that  is 
equivalent  to  an  RMS  uncertainty  in  the  measured  reflectivity  differences 
of  100%.  Clearly,  imaging  and  MTI  systems  cannot  operate  under  these  cir- 
cumstances and  means  must  be  found  to  reduce  this  target  noise. 

It  will  be  shown  that  for  case  where  the  diffuse  signal  dominates  the 
specular  component  the  signal  amplitude  will  obey  Rayliegh  statistics  while 
the  intensity  follows  an  exponential  distribution.  When  a specular  signal 
is  present  with  the  diffuse  component  we  will  show  the  intensity  follows  a 
modified  Rayliegh  or  a Rician  distribution.  In  particular,  it  will  be  sho*..p 
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that  the  specular  component's  signal  to  noise  ratio  will  modify  the  noise 
spectrum  and  control  the  moments  of  the  probability  distribution.  The 
point  that  we  wish  to  emphasize  is  that  the  specular  and  diffuse  components 
are  not  independent  and  the  moments  of  the  diffuse  scintillation  statistics 
depend  on  the  presence  of  the  specular  component.  In  cases  where  the  dif- 
fuse signal  fluctuation  is  high  ard  where  incoherent  addition  of  images  is 
required  it  is  found  that  the  image  improvement  will  depend  on  the  signal 
to  noise  of  the  specular  component. 

The  most  straightforward  way  of  improving  the  imagery  is  to  obtain 
several  images,  each  with  its  independent  noise  characteristics,  and  then 
incoherently  add  these  images.  In  performing  this  addition  the  signal 
increases  as  the  number  of  independent  observations,  N,  while  the  RMS 
fluctuations  increase  as  and  the  image  S/N  improves  as  VnT  Thus,  if 
we  had  an  image  with  a RMS  fluctuation  of  50%  and  we  wish  to  reduce  this 
fluctuation  to  5%  we  would  have  to  obtain  100  dependent  images  to  gain  a 
factor  of  10  in  the  RMS  signal  fluctuation.  If  a system  design  or  mission 
calls  for  a noise  equivalent  reflectivity  difference  of,  say,  .05,  and  if 
we  expect  a diffuse  scintillation  fluctuation  of  50%  we  must  determine  the 
best  way  to  obtain  100  dependent  observations.  It  should  be  emphasized 
that  a large  power  S/N  is  not  required  for  each  independent  observation 
and  this  fact  will  be  of  great  assistance  in  the  system  design  when  avail- 
able transmitter  powers  are  restricted  to  low  values. 

Another  point  to  be  emphasized  is  that  in  order  to  obtain  images  with 
independent  noise  fluctuations  the  target  and  receiver  should  jt  introduce 
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correlations  between  observations.  In  the  sections  that  follow  we  will  de- 
scribe target  scintillation  averaging  techniques  in  which  the  measurements 
are  believed  to  have  correlated  statistics.  It  will  be  necessary  to  per- 
form measurements  with  diffuse  bar  charts  to  demonstrate  that  the  frames 
obtained  are  indeed  independent  and  that  an  improvement  in  the  image  S/N 
increases  as  the  square  root  of  the  number  of  indeoendent  observations. 

When  diffuse  scintillation  reduction  is  necessary  it  will  be  required 
to  combine  images  with  independent  noise  statistics.  The  independence  of 
the  noise  statistics  will  depend  on  the  amplitude  spatial  correlation  at 
the  receiver  aperture.  The  theory  of  this  correlation  is  presented  and  an 
equation  for  the  amplitude  correlation  distance  is  derived. 

We  then  discuss  various  methods  used  to  reduce  diffuse  target  scintil- 
lation which  include:  multiple  observations,  multiple  receiver  and  trans- 
mitter averages,  focal  plane  arrays,  polarization  averaging,  and  frequency 
averaging.  A derivation  of  the  origin  of  the  wavelength  dependence  of  the 
diffuse  scintillation  is  also  presented. 

Finally  we  discuss  methods  whereby  the  speckle  reduction  can  be  ob- 
served in  the  laboratory.  In  particular,  the-  system  noise  equivalent  re- 
flectivity and  angular  resolution  degradation  can  be  observed  by  employing 
diffuse  bar  charts  of  known  surface  characteristics.  From  these  measure- 
ments the  moments  of  the  scintillation  can  be  measured  and  the  decorrelated 
images  can  be  combined  to  observe  the  reduction  in  noise. 

The  author  would  like  to  thank  Or.  David  Fried,  OSC,  Prof.  Nicholas 
George,  C.I.T.,  and  Prof.  Joseph  Goodman,  Stanford  University  for  the  many 
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helpful  discussion  concerning  the  topics  discussed  in  this  report.  Special 
thanks  go  to  Harvey  Buss,  Lincoln  Laboratory, for  his  assistance  with  the 
calculations. 

2.  AMPLITUDE  CORRELATION  DISTANCE 

An  intuitive  description  for  the  origin  of  diffuse  target  scintilla- 
tion may  be  obtained  from  the  following  argument.  Figure  1 shows  a plane 
wave  falling  on  a diffuse  reflector  with  diameter  The  diffuse 

surface  will  introduce  a random  phase  shift  across  the  target  and  in  the 
far  field  a spatially  modulated  field  pattern  will  be  observed.  The  origin 
of  this  spatial  tiv^dulation  is  the  randomly  phase  shifted  Huygen  wavelets 
emanating  from  the  illuminated  surface.  As  the  transmitted  beam  is  scanned 
across  the  diffuse  surface  the  spatial  modulation  of  the  far  field  constantly 
changes  and  the  receiver  will  find  maxima  and  minima  in  a random  fashion. 

The  end  result  is  an  image  with  a highly  granular  appearance  with  a low 
information  signal  to  noise  and  a degraded  angular  resolut’'  i.  The  statis- 
tics of  the  fluctuations  will  be  described  and  it  will  be  shown  that  the 
RMS  fluctuation  can  be  as  large  as  the  signal  mean.  This  can  yield  close 
to  100%  noise  in  cases  where  there  are  small  changes  in  reflectivity. 

In  order  to  perform  diffuse  target  scintillation  reduction  through 
the  incoherent  additions  of  images  with  independent  signal  fluctuations, 
we  need  to  obtain  the  size  of  the  diffuse  target  correlation  length  L^.. 
Consider  the  arrangement  in  Figure  1.  The  diffuse  target  correlation  length 
is  the  distance  required  in  moving  the  receiver  aperture  from  point  to 
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Fig.  1.  A plane  wave  across  a diffuse  surface  creates  a far  field  pattern 
with  a spatial  amplitude  correlation  length  L^.. 
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?2  to  obtain  u signal  at  point  ?2  uncorrelated  with  that  signal  at  . The 
correlation  length  may  be  obtained  by  recognizing  that  it  is  the  separation 
required  to  have  the  amplitude  at  points  1 and  2 uncorrelated. 

The  amplitude  correlation  distance  may  be  obtained  directly  if  one 
treats  the  illuminated  target  as  a spatially  incoherent  source  which  is 
certainly  valid  for  a diffuse  surface  with  a large  number  of  scatterers. 

We  will  show  that  the  phase  correlation  distance  Lj.  is  given  by 


where  R is  the  distance  from  the  target  to  the  position  midway  to  points  1 
and  2.  To  obtain  the  correlation  distance,  we  refer  to  Figure  2 where  the 
source  is  a laser  illuminated  spot  which  we  assume  to  be  a diffuse  surface. 
For  simplicity,  we  imagine  that  the  diffuse  source  is  a collection  of  N in- 
dependent scatterers  and  we  wish  to  obtain  the  amplitude  correlation  at  the 
points  and  P2  in  Figure  2.  is  the  complex  amplitude  at  the  point 
P^  associated  with  the  i^^  scatterer. 

The  total  amplitude  at  P^  is  given  by 


A 


A " 

* 


(1) 


where  I.  is  the  ave-age  intensity  of  the  i^^  scatterer  and  is  its 

range  to  the  point  P^  shown  in  Figu»*e  2. 
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TT-9  (2) 


Fig.  2.  The  geometry  for  a diffuse  source  in  the  plane  (x.y)  and  a receiver 
plane  (X,Y)  at  a range  R. 
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and  the  complex  amplitude  correlation  function  is  given  oy 


- " A >’)  A (2>  * 

1 , j=i  ’ ^ 

- ? ...  N 

■ ,f,  . (1)  . (2)  • + S . (1)  . (2)  * 

' ' "l  "l  i«  N "j 

where  the  bars  denote  an  ensemble  average. 

We  now  employ  our  assumption  of  a diffuse  source  which  implies  that 
the  N scatterers  are  independently  oriented  and  so  large  in  number  that 


iTIT7W 

'i 


0.  m 


and  the  amplitude  correlation  function  simplifies  to 


iWr*  = ? A <’> « * 

i=l  ’ ^ 

= Z I.  exp  [i  k 

i=l  1 1 

If  I (a)  da  is  the  intensity  from  differential  area  da  the  correlation 
function  may  be  written 


5f f I(a)e  ^*^1  ■ da 

A^(Pi)A2  (P2)  =*^3 
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and  we  wil'  now  show  that  this  reduces  to  the  form 


(P^)  Ag  * (Pg)  " J f ^ (x,y)e"'*'‘(px  + qy) 


dx  dy 


where  we  follow  Born  and  Wolf  with 


_ _ -vh  _ -r^ 

P ■ R ’ ^ ' R 


with  R^ , R2»  X ^ , and  shown  in  Figure  2. 

Now  with  (x,y)  the  source  coordinates  and  and  (Xg,^)  as 


shown  we  have 


R,^  ■ (X,  - x)^  + (V,  - y)  ♦ R^ 
^ (X,  - x)^  + {T,  - y)^ 

"l  ' I*  * ?R— ’ 


with  similar  equations  for  Rg^  and  Rg  since  ^^/R,x/R,  and  y/R  are  all  much 
less  than  unity.  The  difference  in  P^  and  Pg  is 

. (Ki'  y/)  - (Xg2  . ^2) 

- Kg  - 2^ 


(X.1  - ^)x  + (Y^  - Yg)y 
R 
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and  substituting  for  p and  q we  have 


Rl  - R2  = 


(X,^  t y/)  - (Xj^  + Yj^) 


- px  + qy 


and  we  write  the  correlation  function  as 


(P^)  Ag  * (P2)  * exp  (iij<)x 


//dxdy  I (x,y)  exp  (-ik(px  + qy) 


where 


x[(x/  t Y,^)  - t 1^^)] 


represents  the  phase  difference  2ir  (0'P-|  - O'Pg)  /X  and  is  usually  neglected 


since 


2TTk(0'P^  - 0’P2)  «ir 


and  we  obtain 


A^  (P^)  A2  * (P2)  =y  y*dxdy  I (x,y) 

X exp  (ik(px  + qy)) 
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'l.!-W|J  HP,  I tjtiUJjSgg'^ 


"OIW 


and  if  the  source  is  a disk  with  radius  p and  I (x,y)  = 
we  obtain 

V 


(P^)  A2  * (P2) 


Vl  <^) 

V 


where 

v=  Vd,  - k>^  + (I,  - 


AR 


where  the  integration  is  the  standard  form  found  in  diffraction  theory. 

is  the  first  order  Bessel  function  of  the  first  kind.  This  result  was 
first  obtained  by  van  Cittert^Hnd  is  the  same  pattern  that  one  would  ob- 
tain if  the  source  was  a diffraction  aperture  where  the  aperture  now  has 
an  amplitude  distribution  which  is  given  by  the  intensity  across  the  source. 

The  spatial  correlation  length  will  be  obtained  when  we  employ  the 
condition  AR  L^,  where  (v)  = 0.  We  have 


2TTkpL 


^ = 1.2tt 


L ^ _i 
‘-C  L 


Targ 


with  L 


Targ 


= 2p. 


' ^ “"^  r 


For  cases  where  the  target  illumination  spot  is  given  by 


'■Targ  = ^ \ 


where  0^^  is  the  transmitter  divergence  we  obtain 


L = 
‘-C 


and  for  a diffraction  limited  transmitter,  0|^  = X/D^,  we  have 


'-C  = ‘>T 


where  Dy  is  the  transmitter  diameter. 

In  sections  5.1  and  5.2  we  will  discuss  the  spatial  and  angular 
averaging  of  images  that  contain  diffuse  target  scintillation  noise.  We 
will  see  that  the  correlation  distance  = Dy  will  be  employed  in  system 
design  to  obtain  images  whose  noise  probability  distributions  are  in- 
dependent. 

3.  WAVELENGTH  DEPENDENCE 

The  origin  of  the  speckle  wavelength  sensitivity  is  the  coherence 
time  of  the  light  source.  If  a light  source  has  a spectral  density  S(a) 
with  an  average  width  of  Ao  (wavenumber,  cm"^)  then  the  temporal  coherence 
time  of  the  light  is  given  by 


where  c is  the  velocity  of  light. 
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We  say  that  the  electromagnetic  wave  consists  of  wavepackets  with  a 
coherence  length  of  cA  t^  Aa~^.  The  physical  interpretation  of  the 
coherence  time  is  the  temporal  analog  to  spatial  coherence  and  has  its 
origin  in  the  amplitude  time  correlation  which  occurs  at  one  point  is  space. 
The  Michel  son  interferometer  spectrometer  measures  this  time  correlation 
(by  introducing  a time  delay  6t  = 6x/c  where  6x  is  the  mirror  retardation) 
and  records  through  the  amplitude  correlation  function 

A (t)  A (t  + 6t)  = A (x/c)  A (x/c  + 6x/c) 

where  the  bars  denote  an  ensemble  average.  The  Fourier  transform  of  the 
interferogram,  A (t)  A*  (t  + 6t),  gives  the  power  spectrum  of  the  field 
S (v)  which  has  a spectral  width  6v  (Hertz). 

Now  consider  a plane  wave  of  spectral  width  Ao  falling  on  a diffuse 
source  which  has  a mean  roughness  height  For  simplicity,  let  us 
consider  the  two  spectral  patterns  which  originate  from  planes  A + B shown 
in  Figure  3.  If  the  distance  I between  planes  A and  B is  much  shorter  than 
the  temporal  coherence  length,  i.e.,  I « Ao~\  then  there  will  be  a 
strong  temporal  phase  correlation  between  the  two  patterns  and  the  speckle 
pattern  from  each  will  be  highly  correlated.  If  however,  i » Aa”\  the 
temporal  phase  correlation  between  the  two  planes  will  be  lost  and  the 
speckle  patterns  will  be  independent.  In  this  case,  the  patterns  from 
plane  A and  B originated  from  wavelengths  which  v/ere  separated  in  time  by 
6t  = %/q  and  this  time  separation  was  large  enough  to  cause  a phase 
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g.  3.  The 
ntribute  to 


decorrelation  in  the  wave. 

We  thus  arrive  at  tha  following  condition.  If  the  RMS  height 
variation  of  a diffuse  source  is  such  that 


2 

the  image  will  be  highly  speckled.  Note  we  have  written  Aa  = AX/X  where 
AX  ts  the  width  of  the  spectrum  in  microns. 

In  the  model  discussed  in  Figure  3 we  have  assumed  that  the  laser 
radiation  was  normal  to  the  planes  A and  B and  that  the  necessary  time 
delay  was  due  to  the  separation  distance  i.  Speckle  pa<-tern  decorrelation 
may  also  be  obtained  from  the  geometry  shown  in  Figure  4 where  the  angle 
of  incidence  is  given  by  the  angle  0 and  the  illuminated  spot  now  takes 
the  shape  of  an  ellipse  instead  of  a circle.  Now  consider  the  simplified 
case  where  the  diffuse  elliptical  reflector  consists  of  one  diffuse  plane 
in  contrast  to  the  dual  plane  source  shown  in  Figure  3.  Due  to  the  fact 
that  sections  of  the  ellipse  have  a different  range  we  can  obtain  independ 
ent  speckle  patterns  from  each  elliptical  cross  secti'^n  which  satisfies 
the  inequality. 


where  6R  is  shown  in  Figure  4 and  is  given  by 


TRANSMITTER/RECEIVER 


V 


L-  <"). 


Fig.  4.  As  in  Figure  3,  a variation  in  range  of  the  scatterers  of  a diffuse 
surface  can  contribute  to  diffuse  target  scintillation  averaging  for 
6R>Aa-1. 


i 


where  r is  the  spot  radius  at  normal  incidence  and  N is  the  number  of  in 

.-1 


dependent  strips  in  the  distance  r (cose)  '.  Thus,  in  order  to  have  N 
independent  speckle  patterns  from  a frequency  source  with  a spectral  width 
Aa  we  must  satisfy. 


N < Aa 


2r 


COS0 


The  methods  to  reduce  the  speckle  contrast  will  be  described  in  the  section 
on  frequency  averaging. 

The  radius  r along  the  x axis  in  Figure  4 is  given  by 


1 


r = ^ R . 0j 


where  R is  the  range  and  0^  is  the  transmitter  divergence  and  for  a 
diffraction  limited  transmitter  we  have 


r = RA/2D.^ 


In  this  case  the  number  of  independent  images  that  will  be  obtained  is 
given  by 


M 1 A., 

N = Aa  • 


D.J.  COS0 


To  obtain  frequency  averaging  during  the  duration  of  one  detector 


dwell  time  tp  the  detector  will  have  to  observe  N independent  phase 


correlations  in  this  interval. 


1 ■} 
< -i 


2 '• 


1 ; 
! 
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That  is. 


f' 

i. 


; 1 


and  an  upper  limit  on  the  number  of  independent  distributions  is  given  by  i 

I 

N = cAo  tg 

1 

Since  for  most  CO2  laser  sources  Ao  is  usually  less  than  a wavenumber  we 
will  not  obtain  frequency  averaging  unless  tp  is  very  long  or  frame  to 
frame  averaging  is  employed.  In  this  case  the  dwell  time  is  replaced  by 
the  image  repeat  time  and  N can  be  considerably  increased. 

4.  DIFFUSE  SCINTILLATION  PROBABILITY  DISTRIBUTIONS 
4.1  No  Specular  Component  Present 

We  consider  the  case  of  a plane  wave  falling  on  a diffuse  target 
which  consists  of  a very  large  number  of  scatterers.  The  complex  amplitude 
is  given  by 


A = A.|  + i A2 
and  the  modulus  by 

A = (A,^  . a/)'/2 

Each  scatterer  will  introduce  a random  shift  to  the  phase  and  in  the  limit 
of  a very  large  number  of  scatterers  we  assume  thut  the  amplitude  will  be 
normally  distributed.  The  joint  probability  density  corresponding  to  the 
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real  and  imaginary  parts  of  the  field  is  given  by 


\ 


P (Ap  A^)  = — ^ 


(2.o2)1/2 


A 2 ^ . 2 

Ai  + A2 


where  a is  the  variance. 


We  wish  to  express  this  probability  density  in  t»rms  of  the  intensity 
and  the  phase.  The  statistics  of  these  new  variables  are  obtained  from 


I = A" 


= A ^ + A ^ 


which  is  the  statistical  variable  which  must  be  considered  when  dealing 
with  a process  such  as  square  law  detection.  The  statistics  of  the  phase 
are  also  important  and  may  be  obtained  from  the  change  in  variables 


(|)  = tan 


A,  = I COS4) 


A,  = sin<|> 


and  the  joint  probability  function  of  I and  (j>  is  then  given  by  the  trans- 
formation 
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P (A^,  A2)  = P cos<j),  Sind)) 

and  since  the  intensity  and  phase  are  independent  we  have 
P (I,  d))  = P (I)  P (d>) 
where  a change  to  polar  coordinates  gives 

P (I)  = -4“  exp  (-I/2a^),  I > 0 
20^^ 


= 0 otherwise 
a ad 

P (<t>)  ^ Yii  ° 

= 0 otherwise 

Thus,  when  no  specular  component  is  present  the  intensity  follows 
an  exponential  probability  distribution  while  the  phase  is  uniformly 
distributed.  When  a signal  is  present  the  intensity  does  not  follow  a 
negative  exponential  distribution  and  will  be  shown  in  a later  section  to 
follow  a Rician  distribution. 

4.2  The  Addition  of  Uncorrelated  Intensity  Distributions 

We  wish  to  obtain  the  new  probability  density  which  is  obtained  from 

* ^ h 

the  addition  of  M speckle  patterns.  We  know  that  the  i*^"  distribution 
follows 
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Consider  first  the  case  where  all  the  original  probability  distribu- 
tions have  the  same  mean  intensity,  i.e.. 


X = Xi  ~ ^2  “ *"  ” 


We  then  obtain 


Pu(I)  = do) 


(1  -M)' 


(M-D!  (xr 


exp  (-I/X) 


In  terms  of  the  normalized  variable  x = I/X  we  obtain 


which  is  shown  in  Figure  5 for  various  values  of  M.  The  variance  of 
this  distribution  is  given  by 


A _ I ^ 


and  the  contrast  improves  as  the  inverse  square  root  of  the  numbers  of 
independent  observations.  Note  we  have  written  (M-1)!  = r(M). 

The  moments  of  this  r-distribution  are  given  by 


■S  =/<"  ‘'phO* 
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f di  I^^'^exp(-I/X) 
< (M-1)!  (A)” 


(X)N  (M+N-1)i 
(M-1)! 


for  A^  = A,  for  i = 1 , . . . , M. 

The  first  and  second  moments  are  given  by 


I(^  = 2a  . M 


= 4o^  M(M+1) 


and  the  variance  is 


~7  4 

■«  - 'm  = 


The  mean  square  noise  contrast  is  given  by 


~T 

r 2/»^  _ ^M  ■ ^M  _ 1 

Cj  (M) 2 M 


I 


M 


and  the  root  mean’ contrast  is 


Ct^'^^(M)  = — 
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which  shows  that  for  M uncorrelated  patterns  the  mean  noise  contrast  im- 
proves as  we  had  earlier  expected. 

We  now  turn  to  the  problem  of  M uncorrelated  speckle  distributions 
which  differ  in  their  mean  intensities,  i.e., 

^1  ^ ^2  ^ * * * * ^ ^[yj 

In  this  case  the  probability  distribution  is 


M 

= I 
i = l 


s- 


M-2 


M 

n (A.-A  ) 

0=1  ^ J 
0^*1 


exp  (-I/A^) 


The  moments  of  this  distribution  are  given  by 


M 

I 

i=l 


M+N-1 

j 

M 


n (x,-xj) 
0=1 
OJ^i 


We  now  take  the  special  case  where  the  mean  intensity  is  given  by 


A.  = A + iAi 
1 0 

This  will  correspond  to  M independent  distributions  which  can  be  arranged 
with  average  intensity  differing  by  AI. 
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I 

1!  i 


I 

I 

( 

! 

! : 


'•  1 


We  then  have 


M 

Pm(I)  = 2 

M i=i 


(Aq  + iAi) 


M-2 


M»-'  n (i-j) 


expj^-I/(A^  + iAI)j 


i7j 


which  is  shown  in  Figures  6-7  for  various  values  of  AA  and  M. 
The  moments  are  given  by 


I 


N 


M 

E 

i=1 


(Aq  + iAi) 


M+N-1 


M 1 M 

(Air*  n (i-j) 
j=l 


It  should  be  pointed  out  that  the  probability  distributions  for  the 
speckle  pattern  are  expected  to  follow  the  distributions  we  have  discussed 
but  this  remains  to  be  observed  experimentally.  To  observe  the  character 
of  these  distributions,  a measurement  of  the  higher  order  moments  should  be 
conducted.  The  higher  order  moments  of  single  patterns  should  first  be  ob- 
tained and  then  these  patterns  should  be  combined  and  the  process  repeated. 
Before  the  individual  patterns  are  combined,  it  should  be  verified  that  they 
are  statistically  uncorrelated.  This  should  be  accomplished  by  calculating 


the  cross  correlation  coefficients  I^. Ij  for  each  position  in  the  image. 
If  a correlation  is  found,  it  should  be  determined  if  this  is  due  to  the 


presence  of  a signal  by  observing  the  spatial  dependence  of  I.I.. 

* \J 

4.3  The  Addition  of  Uncorrelated  Amplitude  Distributions 

So  far  we  riave  been  concerned  with  the  addition  of  speckle  intensity 
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patterns  which  are  obtained  using  energy  detection.  We  now  wish  to  obtain 
the  probability  distributions  for  systems  which  respond  to  the  amplitude. 

We  will  first  derive  the  density  function  for  amplitude  and  then  derive  the 
density  for  M independe  t amplitude  distributions. 

We  wish  to  obtain  the  probability  density  for  the  variable  x = A/o  = 
Since 

p(I)  = — ^ exp  {-I/2a^) 


we  have 


p(x)  = X exp  (“XV2) 


and  the  moments  of  the  normalized  amplitude  are  given  by 


.m 


fr'  dx. 


If  we  let  z = )C/2  we  obtain 


..m 


^ 2'"/2  g-z 


dz 


= 2"’/2  V {^) 


with  the  first  and  second  moments  given  by 


X - /T.  x^=2 


30 


1 


? 


I 

I 
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! 
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r 

i 


t 


which  yields  a square  contrast  factor  of 
= -1  = .27 

which  is  much  less  than  unity  which  was  obtained  for  the  intensity  dis- 
tribution. Note  that  the  above  amplitude  distribution  was  normalized  to 
X = A/a  where  2a  is  the  second  moment  of  A,  i.e.,  this  amplitude  dis- 
tribution has  not  been  normalized  to  its  mean,  but  rather,  to  its  second 
moment. 

The  amplitude  distribution  may  be  normalized  its  mean  by  the 
change  in  variables 


which  yields 


p(y)  = I y exp  (-  J y^) 


with  moments  given  by 


2*^  r 


We  now  have 


y = 1,  <A>  = A 


and 


7-  <A^ 


i 

TI 


31 


r 
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which  yields  the  same  square  contrast  factor  of  = .27. 

We  will  now  show  the  procedure  for  obtaining  the  amplitude  probability 
distribution  for  M independent  observations  when  there  is  no  background 
term  present.  The  technique  is  to  find  the  characteristic  function  for 
one  observation  and  then  to  raise  the  function  to  the  Mth  power.  We  then 
take  the  Fourier  transform  of  the  Mth  power  characteristic  function. 

If  A denotes  the  mean  value  of  the  absolute  value  of  the  amplitude, 
then  the  normalized  variable  y = A/A  yields  the  amplitude  distribution. 

rr  .2, 


p(y)  = J y exp  {- 


where  we  have  explicity  eliminated  the  second  moment  2a  which  was  used  in 
the  intensity  distributions.  The  characteristic  function  corresponding  to 
p(y)  is 

-(+) 

n „ 2* 

y e 


C(w)  = ft'*'  [p(y)] 


I FT<*> 


f r 


= Airiwexp  {-Attw^} 


where  we  have  extended  the  values  of  y to  + ».  When  y is  large  this  does 
not  yield  a large  error.  And  the  Mth  order  characteristic  function  is 
given  by 


5 -S 

I •; 
If  : 


C^(w)  = (4iii)^  w^  expMTrMw^} 
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The  Mth  order  probability  distribution  is  then  given  by 


P„(y)  = ft' 


[.V,] 


I.  ! p'l] 


y » 1 

where  p = 2uiw.  This  Fourier  transform  is  given  in  the  transform  tables 
of  Campbell  and  Foster^^^-  and  is  given  by 


PM(y)  = 


(t)  2 


y » 1 


where  Hu  is  the  Mth  order  Hermite  polynominal  function  and  where  we  have 
M 

extended  the  values  of  y to  + <». 

For  the  more  realistic  case  of  positive  values  for  y we  have 


= FtM  [I  yexpfj/] 

= |-y*y  exp  j^-  I y^  + ZiTiwyj  d 
= 1 [l  -pexpj^j  erfc  j^j] 
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iv 


where  p = 2Triw.  The  FT  of  the  Mth  power  of  c(w)  does  not  exist  in  closed 
form.  Thus,  for  the  case  where  we  wish  to  know  P|y|{y}»  for  all  y > 0,  a 
closed  expression  is  not  obtained.  Although  Pn(|(y)  does  not  exist  in  closed 
form,  we  may  obtain  the  first  few  moments  of  P^Cy)  f<^om  the  moments  of 
p(y)  as  follows.  It  is  shown  in  sec.  4.3  that 


and  since 


^ = 7 + M(M-1)  y ^ 


M _ 2'"  ^ ,M+2x 

^ WIT  ^ 


we  obtain 


^ ^ ^ + M(M-1 ) 


with  a square  contrast  of 

4 


2 _ IT 


-1 


''  M 

4.4  Specular  Component  Present 

We  have  now  seen  that  if  the  number  of  scatterers  is  very  large  and 
uniformly  distributed  the  received  intensity  follows  a Rayleigh  distribution 
in  amplitude  or  a negative  exponential  distribution  in  intensity.  However, 
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these  statistical  distributions  are  not  obtained  if  the  scatterers  are 
nonuni formly  spread  over  the  target  a'ea,  i.e.,  if  the  target  reflection 
is  a function  target  position  or  if  a specular  component  is  present. 

Let  Aj^  represent  the  part  of  the  amplitude  which  is  not  normally 
distributed  and  is  a result  of  the  changing  target  reflectivity.  Then  the 
amplitude  may  be  written  as 

“ 1/2 


A = (A^  + Ap)^  + A2  ^ 


and  we  then  inquire  as  to  the  form  of  the  probability  distribution  of  this 
new  variable. 

The  new  probability  density  has  the  form  of  a Rician  distribution 
and  is  given  by 

P^(A)  - \ exp  1/2  (A^  + A„^)  /o^  J X 
I,  (AA^/e^) 

2 

where  2o  is  the  mean  square  amplitude  of  the  sum  of  the  oscillations 
2 

and  where  (AAj^/a  ) is  a modified  Bessel  function  of  order  zero.  This 
distribution  is  obtained  as  follows.  The  joint  gaussian  distribution  in 
the  presence  of  the  real  signal  is 
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and  since 


2 I r 

d P(A^,A2)  = — 2 — exp  j-  (A^-Aj^cosg, 

]/ 


+ (A2-Aj^sifB  )^  / 2a^ 


2 dA, dA2  I r ? 9 

d P{^^t^2^  = — 2 — exp  I-  A^^  + A2^  - 2A^Aj^cos8 

- 2A2Aj^sinB  + j Za‘ 


A^  = Acosa,  A2  = Asina 


we  change  to  polar  coordinates 


d^P(A,a)  = exp  j " ^ ’ 2AAj^cos(»-B ) + A^^ 


and  the  density  dP(A)  is  obtained  from 
dP(A)  = y da  dP  (A, a) 


^ exp 

a 


(if  + a/)  , , 

-^7-^  I '0  ) 


where 


Iq  (AA^/a^)  = 2h~f^^  ®^P  I 


cosa 
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Note  that  the  original  RMS  noise  was  given  by  o and  we  thus  write 
s = Aj^/a 

as  the  reflected  amplitude  signal  to  noise  ratio.  Let  the  amplitude  be 
measured  in  units  of  o and  define  the  statistical  parameter 

X s A/a 

(4-51 

and  the  density  function  becomes'  ' 

p(x,s)  = X exp  (x^  + s^)/2  j 

This  density  function  is  the  Rician  distribution  and  is  plotted  in  Figure  8 
for  several  amplitude  signal  to  noise  ratios. 

The  first,  second,  and  fourth  moments  are  given  by  Rice^^^  to  be 

X = 1/2  (27r)^/2  g^p  (_52/4)  ^ ^g2)  (^2/^) 

+ s^  (s^/4)j 

? = (2  + s^) 

"T  = (8  + 8s^  + s^), 

A 


P(y) 


The  probability  distribution  in  intensity  is  given  by 


P(I) 


“T  exp 


(I 


I 

r t1/2t  1/2-1 

0 

2 

L 0 

2 

where  is  the  specular  component  in  intensity.  If  r = is  identified 

as  the  intensity  S/N,  and  we  make  all  intensity  measurements  with  respect  to 

2 2 
a we  may  define  y = I/o  and  obtain 


p(y)  = 1/2  exp 


- 1/2  (Y  + r) 


which  is  the  normalized  intensity  distribution  in  the  presence  of  a signal. 
Note  that 


Lim  p(y)  = -U  exp  (-I/2a^) 

which  is  the  negative  expontial  distribution  obtained  in  section  4.1 
p(y)  is  shown  in  Figure  9 for  r = .01,  1,  4,  9,  and  16.  Now,  since 


and 


r = s 
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P(Y) 


-I i.'V  I 


the  first  and  second  intensity  moments  are  given  by 


y = 2 + r 


T 2 

y = 8 + 8r  + r 


with  the  variance 


T T 

y - y 


= 4 (1  + r) 


with  a square  contrast 


cV)= 

' (2  + rf 


which  is  shown  in  Figure  10. 


4.5  The  Addition  of  Independent  Distribution  in  the  Presence  of 
Background 


Let  us  now  consider  the  case  where  we  wish  to  combine  M independent 
noise  patterns  each  described  by  p^(y^r^. ) 


y,-  = 


1 o. 


9 


r,  » I,  / o, 
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The  new  probability  distribution  Pj^(y)  is  given  by 

Ph<V)  = FT"'  [c„(--)] 

where  C^  (w)  is  the  characteristic  function 

M 

C..  = n c.(o)) 

M i=i  1 


c.(oi)  = FT  p^(y^.r^) 


However,  the  Fourier  transform  of  p^(y^.  ,r^)  does  not  exist  in  closed  form 

and  the  distribution  is  difficult  ' ;.,-rf'>s  in  a finite  number  of  terms. 

th 

However,  if  the  M moment  bj^,j  ot  p^(y. ,r^)  is  known  we  can  obtain 

M f'h 

the  moments  a„  from  the  characteristic  function  since  the  moment  of 
m 

P|^(y)  is  given  by 


I 

Ml 


(-1)"’ 


duxii 


C.,(a)) 


o)=0 


= (-1) 


m 


m 


do)' 


m 


M 

n 

j=i 


Cj(“) 


u=0 
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where 


m d 


m 


db) 


c . (w) 
m j'  ' 


w=0 


We 


consider  the  case  where  the  M distributions  p^. (y^ir^)  are  equal.  Then 


(bj)  = C (b)) 


and 


m 


/ , d*”  ^M/  X 

(-')  — i;r-  c (w) 


db) 


m 


b)=0 


The  first  and  second  moments  are 


and 


a^“  = Mb^ 


32^  = Mb2  + M(M-l)b^^ 


and  the  moments  in  the  intensity  are  given  by 


y^  = M(2  + r) 


yJJ^  = M(8  + 8r  + r^)  + M(M-l)  (2+r) 


44 


I 


1 
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with  the  variance 

^ = 4H  (Ur) 

and  the  mean  square  contrast  is 

M(2  + r)^ 

-1  /2 

The  contrast  improves  with  M as  expected. 

4.6  Correlated  Intensity  Distributions 

So  far  we  have  treated  diffuse  scintillation  images  as  having  in- 
dependent statistics  and  we  now  wish  to  investigate  the  more  complicated 
problem  of  patterns  with  a non-zero  correlation.  We  will  treat  the  addition 
of  two  correlated  speckle  patterns  in  the  absence  of  a signal. 

Let  A = + iA^  and  B = + iB^  correspond  to  fields  from  the  two 

patterns.  We  treat  A^ , A2,  B^ , and  as  normally  distributed  variates  with 
equal  weight.  The  fourth  order  joint  probability  density  is  given  by 


1 


P(Ar  A2>  Bp  B^,  p)  = 2 -—2 

(1  - 


exp 


:2,  "1  ' "2  " ‘^l  " “2 


1 r 2 2 2 2 

'a,  + Ao  + b/  + b/ 


2a"(l-p") 


- 2p(A^  B^  + A^  B^)  I 
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j^r  ~~" — <>. 


;(=«<;-.■■’  1 '“-»vir'»i,<_— 


where  p is  the  correlation  coefficient  which  varies  from  one  to  zero 
We  change  variables  to  the  polar  form  of 

A-|  = Acosy*  1^2  ~ ^siny 


(7) 


which  gives 


B-j  = Bcos5 , B2  = Bsin6 


dA^dA2dB^dB2 


J|‘^  dAdBdyd6 


= ABdAdBdyd6 

where  the  Jacobian  is  given  by  (AB)"\ 

The  new  joint  density  is  given  by 


p(A,y;  B,5)  = 


AB 


(2ira^)  (1-p^) 


exp 


* £1 

2o^  ( 


2 (A^  + B^) 

(1-P^) 


- 2pcos(y-6) 

The  joint  distribution  p(A,B)  is  obtained  by  integrating  over  y and  6 from 
0 to  2 :.  Note  that  only  the  phase  difference  y-<5  appears  in  p(A,y;  B,6  ). 
Since  the  reference  phase  is  arbitrary  we  set  y = 0 and  obtain 

2it 

p(A,B)  = JJ  dydBp(A,y;  B,6  ) 
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A^  + B^ 
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a (1-p  ) ( 2a  (1-p  ) 


/ PAB  \ 

lo^d-P^)/ 


where,  as  in  section  4.3,  Iq(x)  is  the  zero  order  Bessel  function  of  the 


with  an  imaginary  argument. 


I„{x)  • J„(ix) 


2-n 


“2"  / 


d(j>e 


XCOS(j) 


The  density  in  intensity  variables 


= A^ 


IS 


Pdl,  I2)  = -J-— — r exp  ] 

’ ^ a^l  - p^)  ( 2a^d-p)) 


I 


1/2  X 1/2 

pll  I2 


° ' 0^(l-p2) 


and  the  density  in  terms  of  the  variable 


(8) 


I = I,  * I2 
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p(I^,  I-I^)  dl^ 


PC)  . / 


= exp 


I ( F(I,a;p) 


where  F(I,a;p)  is  given  by 


F(I,a;p) 


1<?  (1-p)  j 


dl 


1 


The  function  F(I,o;p)  is  not  available  in  closed  form  and  must  be  obtained 
numerically. 

5.  DIFFUSE  TARGET  SCINTILLATION  REDUCTION 

We  will  describe  four  methods  of  obtaining  statistically  independent 
observations  of  diffuse  sources.  These  methods  may  be  summarized  as: 

1.  Spatial  Averaging 

a.  Multiple  observations 

b.  Multiple  receiver  apertures 

c.  Multiple  transmitter  apertures 

2.  Angular  Averaging 

a.  IFOV  overdesign 

b.  Multiple  IFOV  elements 

3.  Polarization  Averaging 

4.  Frequency  Averaging 
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In  this  discussion  that  follows  we  will  assume  the  system  uses  hete- 
rodyne detection  and  that  the  receiver  IFOV  matches  the  transmitter  beam 
divergence.  Since  a heterodyne  receiver  is  diffraction  limited  the  receiver 
aperture  will  be  fixed  once  the  transmitter  divergence  has  been  chosen  from 
a consideration  of  range  and  required  spot  size.  If  direct  detection  was 
employed  one  could  employ  a larger  receiving  aperture  and  thereby  obtain  a 
reduction  in  speckle  directly  through  aperture  averaging.  However,  for 
coherent  systems  this  freedom  does  not  exist  and  one  must  employ  these 
techniques  listed. 

5.1  Spatial  Averaging 
Multiple  Observations 

The  technique  of  multiple  spatial  observations  is  based  on  the  fact 

that  if  the  laser  radar  is  translated  a coherence  length  and  a new  image  is 

(91 

observed  then  this  image  will  have  independent  speckle  statistics.'  ' The 
images  may  then  be  combined  to  reduce  the  scintillation  effects.  If  an 
aircraft  is  moving  over  a scene  with  a speed  V,  then  the  time  required  to 
make  a second  observation  is  given  by  At  = Lj,/V.  For  the  special  case  where 
L-  = Dy  At  is  the  time  required  for  the  aircraft  to  move  the  distance  equal 

V t ) 

to  the  transmitter  diameter.  The  worst  case  occurs  for  slow  moving  air- 

_3 

craft,  3,000  cm/sec,  say,  and  for  Dy  = 10  cm,  we  obtain  At  = 3 x 10  sec, 
which  is  an  order  less  than  a typical  line  scanner  frame  time.  In  this 
case,  multiple  frames  could  be  combined  and  the  major  problems  would  be 
roll  and  vi^'.ational  stability  to  obtain  accurate  image  registration.  For 
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example,  if  ten  frames  were  to  be  combined,  we  would  require  an  angular 
stability  of  one  IFOV  element  over  the  time  it  took  for  the  registration  of 
ten  frames.  To  simulate  these  independent  observations  in  the  laboratory 
one  can  rotate  the  target  with  respect  to  the  laser  axis. 

Another  problem  to  be  addressed  is  the  image  storage  and  eventual 
addition.  Image  storage  tubes  could  be  employed  and  the  storage  and 
additions  using  CCD's  and  CID's  hold  particular  promise. 

It  should  also  be  pointed  out  that  if  many  observations  are  required 
an  increase  in  scanner  speed  may  not  always  be  practical.  For  example, 
with  imaging  systems  that  also  require  high  doppler  resolution,  there  is  a 
minimum  detector  dwell  time  (and,  therefore,  an  upper  limit  on  scan  speed) 
given  by 

t ^ ^ 

where  6v  is  the  required  doppler  resolution. 

Multiple  Transmitter  and  Receiver  Apertures 
A second  technique  using  spatial  averaging  is  to  employ  a multi- 
plicity of  receiver  and/or  transmitter  aperatures  which  are  separated  by 
a correlation  distance  L^,.  In  the  ideal  case  where  the  correlation  dis- 
tance is  a transmitter  diameter  we  would  have  to  maintain  separations  of 
the  order  of  10  cm  for  10.6pm  systems.  This  type  of  separation  is  not 
difficult  to  obtain  and  should  not  cause  major  design  problems.  However, 
multiple  transmitter  and  receiver  optics  must  be  re-aligned  every  time 
the  altitude  is  changed.  This  re-alignment  is  necessary  so  as  to  maintain 
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the  receiver  IFOV  direcu'.y  on  the  laser  illuminated  target.  Techniques  for 
a continual  alignment  with  altitude  have  not  been  developed  and  a feasi- 
bility investigation  should  determine  how  the  S/N  and  angular  resolution 
degrade  in  systems  with  automatic  alignment. 

5.2  Angular  Averaging 
IFOV  Overdesign 

A straightforward  method  for  reducing  speckle  in  a coherent  laser 
line  scanner  is  to  design  the  system  for  more  angular  resolution  than  is 
required  for  a particular  application.  A highly  speckled  pattern  will  be 
obtained  but  since  there  is  more  angular  resolution  available,  one  may 
average  adjacent  high  resolution  elements  and  thereby  obtain  an  angular 
average  with  a corresponding  decrease  in  speckle.  If  photographs  are 
required  the  procedure  is  simple.  Place  the  high  resolution  speckled 
photograph  in  the  out  of  focus  region  of  a camera  and  take  a second, 
slightly  out  of  focus,  photograph.  The  new  photograph  will  have  less 
speckle  and  less  angular  resolution  since  it  was  obtained  with  the  camera 
slightly  out  of  focus.  Of  course,  there  is  no  practical  reason  why  one 
would  want  to  produce  the  second  photograph  because  it  does  not  contain 
any  more  information  than  exists  in  the  first  ohotograph.  However,  for 
systems  that  operate  in  real  time  and  are  to  obtain  automatic  target 
identification  such  a reduction  technique  is  necessary  and  must  be 
accomplished  through  electronic  storage  and  addition. 


The  price  one  pays  for  this  type  of  averaging  is  that  for  constant 
frame  rates  the  IFOV  dwell  time  is  decreased  with  a corresponding  increase 
in  system  bandwidth.  Although  one  increases  the  noise  due  to  the  increase 
in  the  post  detection  bandwidth  one  gains  signal  when  the  adjacent  high 
resolution  elements  are  combined  with  the  result  that  the  final  signal  to 
noise  is  the  same  as  what  would  have  been  obtained  with  an  instrument  de- 
signed for  the  lower  angular  resolution.  Note  that  in  order  to  obtain  a 
smaller  spot  size  one  must  decrease  the  transmitter  divergence.  Similarly 
one  must  use  a smaller  receiver  IFOV  which,  for  diffraction  limited 
coherent  systems,  implies  a larger  receiver  aperture. 

A diversity  of  at  least  two  can  be  obtained  for  imagers  that  use  a 
continuous  raster  scan  across  a diffuse  target.  To  see  this  it  is  only 
necessary  to  compute  the  number  of  independent  illuminated  spots  that  are 
recorded  in  every  dwell  time. 

Let  the  illuminated  spot  diameter  be  given  by  DTay,g-  Then  the 
number  of  independent  speckle  patterns  obtained  in  one  dwell  time  is  given 


where  is  the  distance  scanned  by  the  illumination  spot  in  a dwell  time 
tg.  is  given  by 

0$  = R § tp  = R 0p 
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0 

where  0 is  the  scan  rate  and  0^^  is  the  receiver  IFOV.  For  a laser  trans- 
mitter, the  target  spot  diameter  is 


where  0j 


°Targ  = « ®T 


is  the  transmitter  divergence  and 


and  for  a diffraction  limited  transmitter  and  receiver 


For  a heterodyne  transmitter/ receiver  sharing  the  same  aperture  we  obtain 


Multiple  IFOV  Elements 

In  the  previous  section,  we  found  that  by  using  IFOV  overdesign  and 
combining  adjacent  high  resolution  elements  the  speckled  appearance  of  a 
photograph  could  be  improved  but  its  photon  signal  would  remain  the  same 
Now  consider  a vertical  fan  beam  which  sweeps  the  scene  horizontally.  Let 
this  vertically  formed  spot  be  followed  by  a receiver  with  n^  vertically 
placed  detectors.  With  this  design  there  are  n^  independent  IFOV  elements 
associated  with  just  one  fan  shaped  laser  beam.  Each  IFOV  element  sees  a 
different  part  of  the  laser  illuminated  spot.  If  each  partially  illuminated 
spot  can  be  regarded  as  an  independent  source  with  statistics  uncorrelated 
to  the  other  adjacent  spots  then  each  IFOV  element  will  obtain  an  independ- 
ent set  of  fluctuations.  As  the  aircraft  mounted  sensor  moves  forward  in 
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the  direction  of  the  vertical  beam  each  IFOV  ground  element  may  be  observed 
by  adjacent  detectors  and,  with  appropriate  storage  and  time  delay,  averaging 
may  be  obtained  without  a higher  information  bandwidth. 

The  problems  associated  with  designing  and  operating  a coherent  focal 
plane  array  are  mostly  concerned  with  how  one  matches  one  local  oscillator 
signal  to  many  detectors.  Quadrant  and  five  element  linear  arrays  have  been 
employed  in  coherent  sensors  and  a feasibility  investigation  should  identify 
what  problems  would  be  associated  with  a nine  element  (three  by  three)  array. 

A ten  element  linear  array  with  a matched  fan  also  deserves  careful  considera- 
tion. 

In  applications  where  a multi-element  focal  plane  is  required  it  is 
desirable  to  place  the  detectors  as  close  as  possible.  However,  optical 
cross  talk  should  be  minimized  to  30  or  40  dB  which  will  demand  that  the 
detectors  be  spaced  approximatley  10%  of  the  diffraction  limited  detector 
size.  If  the  spacing  technique  does  not  yield  the  required  cross  talk 
attenuation,  one  will  have  to  consider  an  optically  shielded  plane  array. 

The  coherence  distance  has  been  shown  to  be  given  by 
This  is  the  distance  at  the  aperture  plane  over  which  there  exists  a 
spatial  amplitude  correlation.  Now  consider  a direct  detection  imager 
with  transmitter  and  receiver  aperture  given  by  D.j.  and  Dj^, respectively.  If 
the  transmitter  is  diffraction  limited,  the  target  diameter  is  given  by 

I = R A 

■-Targ 
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and  the  coherence  distance  is  given  by 


D, 


-Targ  "T 

and  the  number  of  independent  coherence  regions  in  the  receiving  aperture 
area  is  given  by 


N is  therefore  the  number  of  independent  speckle  patterns  which  will  bo 
incoherently  combined  at  the  direct  detection  focal  plane.  Thus,  if  the 
receiver  aperture  is  made  larger  than  the  transmitter  aperture,  a spatial 
averaging  will  be  obtained  in  a direct  detection  receiver.  However,  for 
heterodyne  receivers  employing  a single  detector,  this  type  of  aperture 
averaging  is  not  possible. 

The  reason  why  a single  detector  heterodyne  receiver  cannot  be  em- 
ployed in  an  aperture  averaging  system  is  due  to  the  fact  that  heterodyne 
detection  employs  amplitude  reception.  The  size  of  the  effective  aperture 
in  a heterodyne  receiver  is  only  as  large  as  the  region  of  high  amplitude 
correlation  across  the  aperture  plane.  Thus,  if  the  receiver  aperture  is 

made  larger  than  the  coherence  distance  the  receiver  will  have  a signal 

2 

st.’ength  proportional  to  and  not  to  the  physical  size  of  the  receiver. 

Now  this  limitation  of  a heterodyne  receiver  may  be  circumvented 
by  employing  an  aperture  diameter,  D^,  which  is  larger  than  the  coherence 
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distance,  L^,  in  conjunction  with  a focal  plane  array  of  detectors.  Each 
detector  will  be  sized  to  the  diffraction  limit,  i.e., 

IFOV  = -jf- 


and  the  combined  square  array  field  of  view  will  be  matched  to  the  coherence 

region,  i.e., 

0 ^ 
array 


= 4arg/« 

which  is  just  the  subtense  of  the  illuminated  target. 

Each  individual  detector  in  the  focal  plane  array  will  observe  a small 
fraction  of  the  illuminated  spot  as  is  shown  in  Figure  11.  The  target 
observation  diameter  is 

'■'rarg  = 


where  is  the  heterodyne  aperture  which  is  larger  than  the  transmitter 
aperture  diaineter  D^. 

The  observation  spot  will  be  the  source  of  a new  coherence  patch  of 
diameter 

L.  = -7^ — R 

4arg 
= “h 
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but  the  plane  of  this  coherence  patch  will  be  tilted  with  respect  to  the 
aperture  plane  by  an  angle  of  the  order  r^/R  where  r^  is  the  distance  from 

•f  h 

the  center  of  the  target  illumination  area  to  the  i^"  observation  circle 
corresponding  to  the  i^^  element  in  the  focal  plane  array,  r^.  is  shown 
in  Figure  11 . 

Since  all  the  target  subregions  have  a constant  phase  relationship 
to  one  another,  we  must  incoherently  add  the  output  from  each  detector 
to  obtain  a speckle  reduction.  This  would  be  equivalent  to  having  each 
subregion  illuminated  by  lasers  with  no  spatial  correlation.  The  amount 
of  speckle  reduction  is  then  just  the  number  of  detectors  in  the  focal 
plane  array  which  is  given  by 


- I i ^ 

\x/D„  ) -[b^  ) 

which  is  the  same  number  of  independent  speckle  patterns  obtained  in  the 
direct  detection  system.  Note  that  we  have  assumed  an  extremely  close 
spacing  between  each  element  of  the  focal  plane  array.  The  fact  that  we 
have  obtained  the  same  aperture  averaging  as  the  direct  detection  receiver 
should  not  be  surprising  if  one  recognizes  tl^at,  by  the  incoherent  addition 
of  detectors  in  a square  array  (whicn  has  the  same  total  angular  sub- 
tense as  the  direct  detection  system), we  have,  in  effect,  produced  a direct 
detection  receiver  with  the  photon  sensitivity  of  a heterodyne  system. 
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5.3  Polarization  Averaging 

The  laser  target  illumination  is  plane  polarized  whereas  the  reflected 
beam  is  partially  depolarized.  The  partial  target  depolarization  is  depend- 
ent on  the  surface  characteristics  and  very  little  information  exists  on 
target  polarization  properties.  Two  independent  polarized  fields  will  be 
obtained  in  the  far  field  and  these  fields  will  have  independent  diffuse 

target  scintillation  statistics  and  may  thus  be  combined  to  reduce  the 
(9) 

speckle.'  ' The  amount  of  speckle  reduction  depends  on  the  percentage  of 
target  depolarization.  It  is  clear  that  if  a target  converts  50%  of  its 
reflected  photons  into  each  independent  polarization  component  then  a 
diversity  of  two  may  be  obtained  and  the  speckle  noise  may  be  reduced  by 
.707. 

To  obtain  polarization  averaging  two  detector  focal  planes  must  be 
employed.  Each  focal  plane  must  utilize  local  oscillators  with  perpendicular 
polarization  components.  This  appears  to  be  a considerable  hardware  effort 
to  obtain  a maximum  diversity  of  two.  However,  there  are  further  advantages 
to  using  a dual -polarized  receiver.  The  signal  to  noise  of  such  a system 
is  improved  by  a maximum  of  two  since,  in  a single-polarized  receiver,  all 
depolarized  photons  are  not  mixed  with  the  local  oscillator.  A further 
advantage  is  obtained  when  one  recognizes  that  target  identification  and 
MTI  detection  probabilities  can  be  enhanced  since  man-made  targets  will 
often  have  different  polarization  characteristics  than  natural  terrain 
background.  Thus,  given  a mission  which  demands  a probability  of  detection 
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with  a given  false  alarm  level,  the  required  signal  to  noise  can  thereby  be 
reduced. 

It  should  be  pointed  out  that  if  one  wishes  to  investigate  polariza- 
tion effects  on  coherent  imagers  it  is  not  necessary  to  employ  a dual  detec- 
tion system  if  one  employs  a system  with  separate  transmitter  and  receiver. 
With  such  a system,  one  can  obtain  laboratory  data  by  first  obtaining  an 
image  for  a fixed  polarization  and  then  obtain  a new  image  with  the  relative 
receiver/transmission  polarization  rotated  by  90°.  The  two  images  will 
then  have  independent  scintillation  statistics  and  will  yield  the  required 
information  on  target  depolarization.  Rather  than  build  a system  with  a 
dual  polarization  capability,  it  is  first  necessary  to  experimentally  deter- 
mine the  degree  of  depolarization  for  various  targets. 

We  wish  to  obtain  the  probability  of  density  for  a partially 
polarized  diffuse  signal  which  is  obtained  by  a detection  system  which 
responds  independently  to  both  degrees  of  polarization.  Let  the  degree  of 
polarization  by  given  by  P and  the  mean  intensity  in  each  beam  be  given  by 

= 1/2  (1  + P)  I 
I2  = 1/2  (1  - P)  I. 

Since  the  polarization  measurements  are  independent  we  have  for  the 
characteristic  function 

C(z)  = [1  - 0 + P)^z]  [1  - 0 - P)Iz] 
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The  probability  density  function  is  obtained  from  the  Fourier  transform  of 
C(z)  and  we  have 


+ 00 


P(I)  = 
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0+P)  I 2]  [1 
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which  can  be  solved  using  Cauchey's  integral  theorem  giving 


PO)  = 


1 


P I 


jexp/-^^-\  - exp[— ■— 

I \{1+P)T/  \{1-P)I 


Figure  12  shows  P( I)  for  I = 1 and  for  the  polarization  values  of  P = 1.0, 
P = .75,  P = .50. 


5.4  Frequency  Averaging 

In  the  section  on  the  wavelength  dependence  of  the  diffuse  scintilla 
tion  patterns  it  was  shown  that,  due  to  the  temporal  coherence  of  the 
transmitter  signal,  the  pattern  would  be  uncorrelated  provided  that  the 
effective  spectral  width  of  the  source  satisfied 


where  is  the  mean  depth  of  tha  diffuse  target  which  Is  shown  in 
Figure  3.  was  also  shown  to  be  the  effective  size  of  a resolution  ele- 
ment when  the  angle  of  incidence  is  oblique  to  the  target  normal  (see 
(9) 

Figure  4). 

Now  for  an  aircraft  radar  with  an  angle  of  incidence  45°  the  effec- 
tive decorrelation  length  is  given  by  cos  (tt/4)  Ly,  where  Ly  is  the 
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normal  target  size  and  the  required  frequency  separation  is  given  by 


-3  -1 

which  is  5 X 10  cm  for  a target  size  of  one  meter.  This  corresoonds  to 
a frequency  difference  of  150  MHz.  Thus,  if  the  laser  spectrum  has  a width 
of  at  least  150  MHz  and  frequency  bins  are  employed  in  the  receiver,  we  may 
obtain  independent  speckle  patterns  from  each  bin  separated  by  this  frequencv 
difference. 

Consider  the  situation  where  the  transmitter  consists  of  M narrow 
spectral  lilies  of  equal  intensity  but  seoarated  from  one  another  at  fre- 
quency intervals  Av.  We  require  a detector  which  will  be  able  to  discriminate 
between  the  M different  frequencies.  This  could  be  accomplished  bv  using 
the  heterodyne  beat  frequency  of  a local  oscillator  whose  center  frequencv 
was  within  the  power  spectrum  of  the  laser  transmitter. 

The  output  from  each  frequency  bin  would  then  be  incoherently  added 
to  obtain  an  improvement  of  N in  the  speckle  signal  to  noise,  where  M is 
the  number  of  spectral  lines  which  are  effectively  emoloyed.  It  should  be 
pointed  out  that  this  system  only  becomes  oractical  where  very  large  band- 
widths  can  be  tolerated  in  a desigr.  or  where  the  oblique  target  spot  size 
is  much  larger  than  one  meter.  For  example,  in  a ground  maoping  mission 
where  the  oblique  target  size  is  10  meters  then  the  required  bandwidth 
separation  would  only  be  15  MHz  and  considerable  improvements  in  soeckle 
contrast  could  be  obtained  with  a broad  band  laser  source  and  a large 
number  of  frequency  bins. 
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Note  also  that  frequency  averaging  requires  the  application  of  a 
local  oscillator  which  is  independent  of  the  transmitter.  In  some  applica- 


tions this  is  an  added  complexity  which  requires  additional  weight  and 
optical  alignment  requirements. 

6.  Diffuse  Target  Scintillation  Measurements 

A summary  of  Rayleigh  and  Rician  statistics  has  been  given  and  we 
wish  to  emphasize  here  that  the  expected  statistics  must  be  verified  ex- 
perimentally by  measuring  the  higher  moments  of  the  fluctuation.  It  should 
also  be  demonstrated  that  the  correlation  distance  satisfies  Lj,  = 
and  that  the  moments  of  the  fluctuations  are  independent  of  provided 
that  there  exist  many  independent  scatters  in  the  target  area. 

The  moments  of  the  fluctuations  should  be  measured  for  a variety  of 
sources  to  see  if  there  is  a difference  in  the  far  field  amplitude  pattern 
for  different  surfaces.  The  correlation  distances  should  also  be  measured 
for  various  diffuse  sources  and  the  inverse  dependence  of  the  correlation 
distance  on  the  target  spot  size  should  be  verified,  that  is,  Lp  = 


I* 

For  each  diffuse  source  we  will  also  need  to  know  the  ratio  of  diffuse 

do' 

to  specular  reflectivity  and  if  the  reflectivity  varies  across  the  surface.'  ' 
From  this  preliminary  investigation  of  diffuse  materials  a few  must  be 
selected  to  construct  the  appropriate  bar  charts  which  we  will  now  discuss. 

The  quantitative  evaluation  of  coherent  images  will  be  obtained  from 
a diffuse  bar  chart  whose  surface  requirements  are  to  be  determined 
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experimentally.  The  type  of  measurements  that  must  be  made  with  the  diffus 
bar  chart  area  fall  into  the  categories  of  the  system  NEAp  and  image 
resolution. 

The  coherent  system  noise  equivalent  .eflectivity  degradation  must 
be  measured.  This  may  be  accomplished  by  using  bar  chart  backgrounds  which 
are  diffuse  and  differ  in  average  reflectivity.  The  bar  chart  must  be 
separated  into  two  parts:  target  and  variable  reflectivity  background. 

As  the  reflectivities  of  the  bar  targets  and  background  are  brought  closer 
together  the  percentage  of  speckle  noise  will  be  obtained  directly  from 
the  NEAp.  Improvements  in  the  NEAp  should  be  observed  in  the  speckle  re- 
duction experiments  which  will  be  discussed  in  a later  section. 

The  coherent  sensor  angular  resolution  degrades  with  increasing 
diffuse  target  scintillation.  While  the  NEAp  experiments  are  being  per- 
formed the  angular  resolution  of  the  system  may  be  obtained  from  the 
diffuse  bar  chart  which  is  constructed  of  diffuse  bars  of  varying  spatial 
separation.  The  final  result  of  the  experiment  will  be  the  graphical  pre- 
sentation of  angular  resolutions  dependence  on  NEAp.  Once  this  type  of 
relation  is  obtained  for  different  diffuse  sources,  one  may  then  determine 
the  required  quantitative  speckle  reduction  for  a system  to  yield  a given 
angular  resolution  and  NEAp. 

An  investigation  should  be  conducted  to  determine  if  the  diffuse 
surface  fluctuations  distort  the  edges  of  man-made  objects.  For  example, 
it  is  not  known  if  a diffuse  background  will  create  granular  effects  along 
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a straight  edge.  This  edge  distortion  appears  qualitatively  in  coherent 
images  but  its  origin  is  poorly  understood.  This  effect  can  be  seen  when 
imaging  a plane  on  a runway.  Often  the  edge  of  the  winds  have  a rough  tex- 
ture and  it  is  not  known  if  this  noise  is  due  to  edge  diffraction,  glints, 
or  diffuse  fluctuations.  It  would  be  a straightforward  procedure  to  place 
straight  edged  objects  above  a diffuse  background  and  determine  the  degree 
of  the  diffuse  scintillation  contribution. 

After  the  contrast  and  angular  resolution  measurements  have  been 
performed,  it  is  then  necessary  to  demonstrate  that  the  speckle  may  be 
reduced  by  a variety  of  techniques.  In  particular,  it  should  be  demonstra- 
ted that  by  moving  the  receiver  a correlation  distance  and  obtaining  a new 
image  pattern  one  obtains  a reduction  in  speckle  contrast  with  an  improve- 
ment in  angular  resolution.  The  improvement  in  contrast  and  angular  resolu- 
tion should  be  graphed  as  a function  of  the  number  of  observations.  If  the 
contrast  does  not  improve  as  the  square  root  of  the  observation  number, 
then  the  observations  are  partially  correlated  and  the  receiver  correlation 
distance  must  be  increased  and  the  experiment  should  be  repeated.  In  order 
to  perform  the  measurements  which  have  been  discussed  it  is  necessary  to 
employ  a laser  line  scanner  with  characteristics  suitable  for  a scintilla- 
tion reduction  investigation. 

The  instantaneous  field  of  view  for  the  coherent  sensor  should  be 
between  .05  mrad  and  .5  mrad.  The  .5  mrad  IFOV  will  be  subject  to  less 
atmospheric  scintillations  in  experiments  where  the  transmit/receiver  path 
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is  greater  than  a few  kilometers.  Since  the  aperture  of  the  receiver  is 
diffraction  limited  the  receiver  aperture  is  dependent  on  the  chosen 
instantaneous  field  of  view.  The  receiver  aperture  diameter  is  then  given 
by 

where  the  IFOV  is  given  by  the  e~^  point  of  the  diffraction  pattern  at  the 
focal  plane.  For  an  IFOV  of  .1  mrad  we  have  a required  aperture  diameter 
of  6.4  cm.  If  an  investigation  using  a focal  plane  array  is  desired,  it 
would  be  preferable  to  employ  as  small  an  IFOV  as  possible  without  being 
distrubed  at  atmospheric  turbulence.  A practical  upper  limit  for  a receiver 
aperture  diameter  is  15  cm  for  ranges  of  a few  kilometers.  If  long  range 
measurements  are  desired.it  would  be  advisable  to  keep  the  maximum  aperture 
to  no  more  than  15  cm  or  20  cm. 

The  transmitter  aperture  should  be  separate  from  the  receiver  so 
spatial  correlation  measurements  can  be  conducted.  Moreover,  the  distance 
between  the  transmitter  and  receiver  aperture  should  not  be  fixed  and  a 
maximum  separation  of  30  cm  should  be  possible.  Note  that  a boresight 
alignment  fixture  will  be  required  since  a separated  receiver  and  trans- 
mitter must  be  tilted  with  respect  to  another  if  they  are  to  observe  the 
same  target. 

To  be  able  to  change  the  target  illumination  diameter  it  will  be 
necessary  to  have  a transmitter  with  a variable  beam  divergence.  This 
could  be  accomplished  by  having  a variety  of  beam  expanders  available. 
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These  beam  expanders  should  mount  to  the  instrument  frame  in  a fashion  that 
will  allow  a continual  boresight  adjustment  with  the  receiver.  The  beam 
divergence  capability  should  vary  from  .1  IFOV  to  10  IFOV  (IFOV  refers  to 
the  receiver's  instantaneous  field  of  view). 

Experiments  involving  specular  targets  wil'  require  that  the  receiver 
be  capable  of  receiving  signals  at  least  10^  times  the  quantum  noise  limit. 
Moreover,  the  system  response  over  this  dynamic  range  should  be  linear  to 
within  10%.  If  experiments  are  carried  out  at  close  ranges  where  the 
average  signal  is  high,  then  the  system  should  have  neutral  density  filters 
available.  The  neutral  density  filters  should  be  capable  of  an  attenuation 
of  10"^  with  grades  of  10"^. 

The  test  system  will  be  required  to  perform  imaging  as  well  as 
moving  target  indication.  The  MTI  minimum  resolvable  velocity  difference 
is  given  by 


6 V = 


fpA 

2 


where  f^  is  the  doppler  bandwidth.  For6v  = 3m/sec,  we  obtain  a doppler 
bandwidth  of  .6  x 10®  Hz  which  gives  a detector  dwell  time  of 

1 _g 

t.  = = .8  X 10~  sec. 

d 2fp 
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